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We present a measurement of the direct CP asymmetry in the decay B° — > K + -k~ using a 
data sample of 227 million T(4S) — > BB decays collected with the BABAR detector at the PEP-II 
asymmetric-energy e + e~ collider at SLAC. We observe a total signal yield of n K - 7T + + n K +^- = 
1606 ± 51 decays and measure the asymmetry (n K - 7T + — n K + 7T -) / (n K - 7I + + n K + 7T -) — —0.133 ± 
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0.030 (stat) ± 0.009 (syst). The probability of observing such an asymmetry in the absence of direct 
CP violation is 1.3 x 10 _J , corresponding to 4.2 standard deviations. 

PACS numbers: 13.25.Hw, 11.30.Er 12.15.Hh 



CP violation has been established in processes involv- 
ing B°-B° oscillations through measurements of the time 
dependence of neutral- B-meson decays to final states 
that include charmonium 0, Direct CP violation, 
a phenomenon that does not involve particle- antiparticle 
oscillations, has been observed in K® decays Q, where 
the effect is a few parts per million. In contrast, a large 
effect is expected in the B-meson system if CP violation 
arises from the Kobayashi-Maskawa quark-mixing mech- 
anism 0,0- 

The Belle Collaboration has reported evidence of di- 
rect CP violation in the decay B° — > 7r + 7r~ at the level 
of 3.2a 6], though this is not confirmed by our measure- 
ment based on a significantly larger data set 0. In this 
Letter we report a measurement of direct CP violation 
in the decay B° — > K + ir~ [j| at the level of 4.2c using a 
sample of 227 million BB pairs collected with the BABAR 
detector at the SLAC PEP-II e + e~ asymmetric-energy 
storage ring. 

Direct CP violation is observable as an asymmetry in 
yields between a decay and its CP conjugate when at 
least two contributing amplitudes carry different weak 
and strong phases. In the standard model, the de- 
cay B — > K + ir~ occurs through two different mecha- 
nisms ( "penguin" and "tree" ) , which carry different weak 
phases and, in general, different strong phases. The di- 
rect CP- violating asymmetry || is defined as 

AffTr = ; , (1J 

n K --K+ + n K+7T - 

where n K -^+ and n K + 7l - are the measured yields for 
the two final states. The charge of the kaon identifies 
the flavor of the decaying B meson (B° — > K + ir~ , B° — > 
K~it + , neglecting second-order weak transitions). The 
Belle collaboration recently reported a measurement of 
Ak-k = -0.088 ± 0.035 ± 0.013 using a data sample of 
152 million BB pairs |l0j . which agrees with our previous 
result 11], and with a less- precise measurement from the 
CLEO collaboration [l^. 

The BABAR detector is described in detail else- 
where [l^. The primary components used in this anal- 
ysis are a charged-particle tracking system consisting of 
a five-layer silicon vertex tracker (SVT) and a 40-layer 
drift chamber (DCH) surrounded by a 1.5-T solenoidal 
magnet, an electromagnetic calorimeter (EMC) compris- 
ing 6580 CsI(Tl) crystals, and a detector of internally 
reflected Cherenkov light (DIRC), providing K-n sepa- 
ration over the range of laboratory momentum relevant 
for this analysis (Fig.E). 

We reconstruct two-body neutral- B decays from pairs 
of oppositely-charged tracks located within the geomet- 



ric acceptance of the DIRC and originating from a com- 
mon decay point near the interaction region. We re- 
quire that each track have an associated Cherenkov-angle 
(8 C ) measured with at least five signal photons detected 
in the DIRC. The value of 8 C must agree within 4<r 
with either the pion or kaon particle hypothesis. The 
last requirement efficiently removes events containing 
high- momentum protons. Electrons are explicitly re- 
moved based on energy-loss measurements in the SVT 
and DCH, and on a comparison of the track momen- 
tum and associated energy deposited in the EMC. We 
use the 8 C measurement to separate kaons and pions in a 
maximum-likelihood fit that determines signal and back- 
ground yields corresponding to the four distinguishable 
final states (ir + ir~, K + tt~, K~it + , K + K~). 

Signal decays are identified using two kinematic vari- 
ables: (1) the difference AE between the energy of the 
B candidate in the e + e~ center-of-mass (CM) frame 
and \/s/2 and (2) the beam-energy substituted mass 
toes = \J (s/2 + pi • p B ) 2 /Ef - p|. Here, y/s is the to- 
tal CM energy, and the B momentum pe and the four- 
momentum of the e + e~ initial state (£Jj,pi) are defined 
in the laboratory frame. For signal decays, AE and Toes 
are distributed according to Gaussian distributions with 
resolutions of 27MeV and 2.6MeV/c 2 , respectively. The 
distribution of Toes peaks near the B mass for all four 
particle combinations. To simplify the likelihood defini- 
tion, we reconstruct the kinematics of the B candidate 
using the pion mass for both tracks. With this choice, 
B° — ► 7r + 7r~ decays peak near AE = 0. For B decays 
with one or two kaons in the final state, the AE peak 
position is shifted and parameterized as a function of the 
kaon momentum in the laboratory frame. The average 
shifts with respect to zero are —45 MeV and —91 MeV, 
respectively. We require 5.20 < toes < 5.29GeV/c 2 and 
|Ai?| < 150 MeV. The large sideband region in toes is 
used to determine background-shape parameters, while 
the wide range in AE allows us to separate B decays to 
all four final states in the same fit. 

We have studied potential backgrounds from higher- 
multiplicity B decays and find them to be negligible in 
the selected AE region. The dominant source of back- 
ground is the process e + e~ — > qq (q = u,d,s,c), which 
produces a distinctive jet-like topology. In the CM frame 
we define the angle 8s between the sphericity axis Hi of 
the B candidate and the sphericity axis of the remaining 
particles in the event. For background events, |cos#s| 
peaks sharply near unity, while it is nearly flat for sig- 
nal decays. We require |cos6>s| < 0.8, which removes 
approximately 80% of this background. 

The selected sample contains 68030 events and is com- 
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FIG. 1: (a) Cosine of the polar angle (8) as a function of laboratory momentum for kaons and pions in simulated B° — > K + tt~ 
decays at BABAR. At a symmetric e + e~ collider operating at the T(4S) resonance, particles from two-body B decays are nearly 
monoenergetic with p ~ 2.6 GeV/c. The boost at PEP-II results in an approximately uniform distribution of laboratory momenta 
between 1.5 and 4.5 GeV/c, and induces a correlation between momentum and polar angle, (b) The measured Cherenkov angle 
for pions (upper band) and kaons (lower band) from D* + — > D°tv + , D° — > K~-k + decays reconstructed in data. The curves 
show the expected angle 9 C as a function of laboratory momentum, for the K and n mass hypothesis, (c) The average difference 
between the expected value of 9 C for kaons and pions, divided by the uncertainty, as a function of momentum. 



posed of two-body B decays (signal) and combinations of 
real kaons and pions produced in qq events (background) . 
We use an unbinned, extended maximum-likelihood fit 
to extract yields, the signal asymmetry {Akk), and the 
background asymmetry {A Kn ). The fit uses mres, AE, 
6> c , and the Fisher discriminant T described in Ref. PJJ to 
distinguish signal and background components for each 
of the four 7r + 7r~, K + ir~, K~ir + , and K + K~ combina- 
tions. The likelihood for event j is obtained by summing 
the product of the event yield and probability "Pi over 
the signal and background hypotheses i. The total like- 
lihood for the sample is 



£ = exp I - ^2 n i El 



(2) 



The probabilities Vi are evaluated as the product 
of the probability density functions (PDFs) with pa- 
rameters d*j, for each of the independent variables 

xj = {tues, AE,T, ®t i ^cT}' wnere ®t an< ^ ®c are the 
Cherenkov angles for the positively- and negatively- 
charged tracks, respectively. We have verified that there 
are no significant correlations between the Xj. For both 
signal and background, the K^ir^ yields are parameter- 
ized as iijc± t t = tik-k (1 T Ak-k) /2, and we fit directly 
for the total yield tlkk and the asymmetry Ak-k ■ 

The 8 C PDFs are obtained from a sample of approxi- 
mately 430000 D*+ -> D°ir+ (D° -> K~tt + ) decays re- 
constructed in data, where K T /n^ 1 tracks are identified 
through the charge correlation with the tt^ from the D*^ 
decay. Figure H^b) shows the measured values of the 
Cherenkov angle as a function of laboratory momentum 
for tracks from the D* sample, and the expected values 
for kaons (8^) and pions (6™ ). Figure [Tfc) shows the av- 



erage K-tt separation, defined as — 6* | /a$ c , where 
ag c is the average uncertainty for kaon and pion tracks 
for a given momentum. The PDFs are constructed sep- 
arately for K + , K~ , 7r + , and 7r~ tracks as a function of 
momentum and polar angle using the measured and ex- 
pected values of 9 C , and its uncertainty. We use the same 
PDFs for signal and background events. 

A total of 21 parameters are varied in the fit. Signal 
and background yields and Kit asymmetries are deter- 
mined simultaneously with the parameters of the signal 
PDFs for rriES and AE, as well as the background PDF 
parameters for toeSi AE 1 , and T . The parameters de- 
scribing the signal T distribution are fixed to the values 
obtained from a large sample of simulated events, and 
the parameters of the C PDFs are fixed to the values ob- 
tained from the D* study. The analysis was performed 
with the value of Aktt hidden until the event selection 
and PDF definitions were finalized. 

The fitted signal yields are uk-k = 1606 ± 51, n„ = 
467±33, and ukk — 3± 12, which are all consistent with 
our previously published measurements of the flavor- 
averaged branching fractions in these decay modes [llj| . 
The direct CP- violating asymmetry is 

Ak-k = -0.133 ± 0.030 (stat) ± 0.009 (syst), (3) 

and the background asymmetry is A\ T , — 0.001 ± 0.008. 
This result is consistent with, and supersedes, our pre- 
vious measurement ^l|. The correlations of Akk with 
A h KlT and ukk are —8% and +2%, respectively. Corre- 
lations with the remaining free parameters are all 1% or 
less. 

The dominant source of systematic error is the poten- 
tial difference between kaons and pions in the dependence 
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FIG. 2: Distributions of AE in data (points with error bars) 
and the PDFs (curves) used in the maximum likelihood fit 
for K + %~ (solid circles and solid curve) and K~tv + (open 
circles and dashed curve). The data are weig hted using the 
background-subtraction technique of Ref. [13 (see text). 



FIG. 3: (a) Distribution of tties enhanced in K + n~ (solid 
histogram) and K~-k + (dashed histogram), (b) Asymmetry 
Ak-tt calculated for ranges of mES- The asymmetry in the 
highest mES bin is somewhat diluted by the presence of back- 
ground. 



of track reconstruction and particle identification on the 
charge of the particle. To estimate this systematic uncer- 
tainty, we use the statistical uncertainty (0.008) on the 
measurement of A\^ as a conservative systematic error 
on Aktt- This background is due to combinations of real 
kaons and pions in the same momentum and polar-angle 
range as the signal tracks, and should have similar sensi- 
tivity to a potential bias. We have also investigated po- 
tential differences in efficiencies for track reconstruction, 
and for the requirement of a minimum number of signal 
photons detected in the DIRC. Using the large sample 
of kaons and pions from the D* study, we confirm that 
the efficiency asymmetries between K + / K~ and ir + /tt~ 
are consistent with zero within the small error of the 
measurements (0.002). Doubly-Cabibbo-suppressed D° 
decays (D° — ► K + ir~) would produce a bias in the 6c 
PDFs derived from the D* sample, but are a negligible 
effect given the current size of the data set. 

We confirm that we are sensitive to a nonzero value 
of Ak-tt by performing fits on samples of Monte-Carlo- 
simulated signal events, and background events gener- 
ated directly from the PDF shapes. With a generated 
asymmetry of —10%, the average fitted value in the en- 
semble of events is Ak-k = -0.102 ± 0.002. Although the 
result is consistent with the generated value, we take the 
sum in quadrature of the error and the difference with re- 
spect to the generated value as a systematic uncertainty 
(0.003). The systematic errors from uncertainties in the 
distribution of T for signal events (0.001) and from the 
parameters describing the 8 C PDFs (0.001) are negligi- 



ble. The total systematic error (0.009) is calculated as 
the sum in quadrature of the individual uncertainties. 

Figure |21 shows background-subtracted distributions of 
AE for signal K + tt~ and K~ir + decays. The subtraction 
is performed using the technique described in Ref. |l5l |. 
where each event is given a statistical weight that de- 
pends on the PDFs and covariance matrix from a fit ex- 
cluding the variable being plotted. The resulting distri- 
bution is normalized to the signal yield and its shape can 
be compared with the PDF we use in the full fit. We see 
no evidence of an enhancement near AE = 0, which could 
arise from significant contamination of B° — > tt + it~ de- 
cays due to imperfect parameterizations of the 8 C PDFs. 

As a further consistency check on the fit result, in 
Fig. we show distributions of m^s for samples en- 
hanced in signal Kir decays using probability ratios based 
on the PDFs for AE, J 7 , and 9 C . The efficiency of the se- 
lection is approximately 80% for signal Kn decays, while 
the contamination from B° — > ir + ir~ is less than 2%. 
Fig. Ofb) shows the resulting distribution of Ak-k as a 
function of mES- 

A number of consistency checks are performed to vali- 
date the result. We generate and fit a large set of pseudo- 
experiments, where the variables Xj for each event are 
generated randomly from the PDFs, and confirm that 
the value of Ak-k is intrinsically unbiased. To check for 
a potential effect from K—tt misidentification, we fit the 
subsample of events (less than half) where both tracks 
have laboratory momentum less than 3.5 GeV/c. The K— 
7T separation for all tracks in this sample is greater than 
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TABLE I: The signal yield uk-k, and signal (Ak-k) and back- 
ground (-4xtt) asymmetries measured in different data-taking 
periods. The number of BB pairs N B g (in millions) for each 
data set is also given. 



Sample 






Akk 


A h K7r 


1999-2001 


21.1 


142 ± 15 


-0.240 ± 0.102 


0.006 ± 0.026 


2002 


66.4 


479 ± 27 


-0.102 ±0.055 


-0.008 ±0.015 


2003 


34.1 


241 ± 19 


-0.109 ±0.079 


0.007 ±0.021 


2004 


104.9 


743 ± 33 


-0.142 ± 0.044 


0.004 ±0.012 



3cr, and we find Ak-k = -0.151 ± 0.047. We perform 
a B —B mixing analysis on the full two-body sample 
using £>-flavor identification and decay-time information 
as described in Ref. From this fit, we simulta- 

neously determine the B lifetime tb = 1-60 ± 0.04 ps, 
B°-B° mixing frequency Arrid = 0.523 ± 0.028 ps -1 , and 
Ak-k — —0.126 ± 0.029, where the errors are statistical 
only. The values of tb and Ama are consistent with the 
world averages |16( , and Akk is consistent with the nom- 
inal fit, demonstrating that the signal events have the 
expected time evolution. Finally, we divide the full sam- 
ple into the approximate period in which the data were 
recorded (Table |]}. We find Akk < and a background 
asymmetry consistent with zero in each data set. 

The statistical significance of the measurement (4.3ct) 
is computed by taking the square root of the change 
in 21n£ when Akk is fixed to zero. If we include 
the systematic error by summing in quadrature with 
the statistical uncertainty, the significance is 4.2cr, and 
the probability of obtaining a negative asymmetry of 
this magnitude or larger in the absence of CP violation 
is 1.3 x 10 ~ 5 . We conclude that the measurement of 
Akk = -0.133 ±0.030 (stat) ±0.009 (syst) reported here 
represents compelling evidence for direct CP violation in 
the i?°-meson system. 
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